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ABSTRACT 



A brief introduction to the gauge invariant classical and quantum theory of cosmological 
perturbations is given. The formalism is applied to inflationary Universe models and yields a 
consistent and unified description of the generation and evolution of fluctuations. A general 
formula for the amplitude of cosmological perturbations in inflationary cosmology is derived. 
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1. Introduction : 

According to the cosmological principle, the Universe should be homogeneous on large 
scales. The isotropy of the cosmic microwave background temperature to an accuracy of 
better than 10~^ is a powerful confirmation of this principle. As a point of further support, 
the most recent large-scale redshift surveys^-* indicate a convergence to homogeneity also in 
the distribution of light. 

However, on smaller scales inhomogeneities exist: galaxies, cluster of galaxies, voids and 
superclusters. The isotropy of the microwave background on smaller scales is an imprint of 
the homogeneity of the matter distribution at the time of recombination. Hence, it is rather 
natural to work under the hypothesis that the present structure of the Universe originates 
from the growth of initially small cosmological perturbations. 

At first glance, the theory of linear cosmological perturbations appears straightforward. 
Given a Priedmann- Robertson- Walker (FRW) background model ^gj^J , Tj^l^^ and small per- 
turbations {Sg^i, , (5T^iy) of metric and energy-momentum tensor, we linearize the Einstein 
equations 

Gfjiv ^ SnGTfj,,, (1) 
about the background solution to obtain 

SG^^ = SnGST^„ . (2) 

The goal of the analysis of these equations is to find the time dependence of the fractional 
density contrast Se/e. 

Linear cosmological perturbation theory was first developed by Lifshitz in 1946, but prior 
to 1980 there was missing motivation for any in depth study, the reason being that there was 
no causal theory for the origin of fluctuations and hence no reason to study perturbations 
except on length scales smaller than the Hubble radius where Newtonian theory is adequate. 

With the advent of inflationary Universe models, the situation changed drastically. As 
shown in Fig. 1, provided that the period of inflation is sufficiently long, all scales of 
cosmological interest originate inside the Hubble radius during the de Sitter phase. Since 
there is Hawking radiation in the de Sitter space^^ with temperature Th ~ H, where H is 
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Figure 1: Evolution of scales in the inflationary Universe. The comoving wave- 
length A of a perturbation is constant in comoving coordinates Xc- The Hubble ra- 
dius increases after inflation {t > tji), but decreases exponentially during inflation 
(t < tn). 

the expansion rate, fluctuations are produced. These inhomogeneities evolve on scales much 
larger than the Hubble radius. Hence, a general relativistic analysis is required. 

The Hawking radiation argument for the origin of perturbations given above is too 
naive^) . The correct analysis uses the familiar quantum particle production effects for quan- 
tum matter fields in an expanding background, as applied to the scalar fields which drive 
inflation. If quantum fluctuations provide the seed perturbations for structure in the Uni- 
verse, then a quantum theory of cosmological perturbations is required. 

Hence, the key issues within the theory of cosmological perturbations are 

- to understand the growth of inhomogeneities on scales larger than the Hubble radius, 

- to develop a quantum theory of cosmological perturbations, 

- to explain the quantum to classical transition for fluctuations. 

In this lecture, we will develop the classical and quantum theory of cosmological per- 
turbations, based on a recent comprehensive review article^'*. The issue of the quantum to 
classical transition will not be addressed (see e.g. Ref. 5 for literature on this topic). First, 
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we will demonstrate that the analysis of the perturbation equation (2) is not trivial: there 
are gauge ambiguities, and the best way to address this problem is to adopt an explicitly 
gauge invariant formalism. 

2. Issues of Gauge : 

In a general context, the gauge ambiguity can be described in two ways. In the passive 
view, we are given a space-time manifold M, a physical quantity Q defined on M, and some 
corresponding coordinate function ^'^^Q{x) (in the above example, (°)(3(x,t) = (t/to)^/^)- Let 
us now choose two sets of coordinates x and x on M. For the first choice, the perturbation 
5Q{p) of (5 at a point p e M is defined as 

5Q{p) = Q{p) - ^'^Qi<p)) , (3) 

whereas for the second choice 

SQ{p) = Q{p) Q{x{p)) . (4) 

For small coordinate changes, the transformation 6Q{p) — > SQ{p) is called a gauge transfor- 
mation. 

In the active view, we are given two manifolds, the space-time manifold M and an un- 
perturbed background manifold N with a fixed coordinate choice. To each physical quantity 
Q on M there is a corresponding function ^^'^Q{x) on N. Any coordinate choice on M corre- 
sponds to some map from N to M (see Fig. 2), and hence to a different definition of SQ{p) 
(see (3) and (4)). 

There are two approaches to the gauge problem. One is to fix the gauge, the other is to 
work in terms of gauge invariant variables. We will now argue that the use of gauge invariant 
variables has many advantages. 

Gravity is not the only theory with gauge ambiguities: electromagnetism is another 
important example. In electromagnetism we can either work in terms of the gauge dependent 
potential or in terms of the gauge invariant field strength tensor F^i,. When using A^, 
the homogeneous Maxwell equations are automatically satisfied, and only the inhomogeneous 
ones need to be solved explicitly. Thus, working in terms of ^4^ makes the analysis easier in 



4 



Figure 2: The active view of a coordinate transformation: the two mappings V 
and V from the background manifold N to the physical manifold M give rise to 
two different coordinatizations of M and hence to differing definitions of perturbed 
quantities. 

the sense that less equations must be solved. The disadvantage is that the variables have no 
direct physical meaning, and that gauge artifacts like the Dirac string may appear. 

In gravitational perturbation theory, however, no simplification of the equations is 
achieved by using gauge dependent variables. Rather, there are more equations and the 
analysis is more difficult. The interpretational problems remain. Hence, there is strong 
motivation to adopt the gauge invariant formalism. 

There is an additional reason for favoring the gauge invariant approach over working in 
the usual gauge- synchronous gauge. In synchronous gauge there is a residual gauge free- 
dom which leads to unphysical modes. Although it is in principle possible to subtract these 
modes, in practice there are formidable difficulties, especially when working with approxi- 
mate solutions. 

Early attempts to develop a gauge invariant theory of cosmological perturbations go 
back some time^\ The first completely gauge invariant analysis was achieved by Gerlach 
and Sengupta^-* and Bardeen^^ . This lagrangean approach was further developed and clarified 
in several papers^). More recently, an alternative Eulerian (or covariant) analysis has been 
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developed in Refs. 10 and 11 and in many subsequent papers^^^. The equations of Ref. 8 
were rederived in Ref. 13 using the Arnowitt-Deser-Misner approach. For a recent review 
of the classical and quantum theory of cosmological perturbations, the reader is referred to 
Ref. 4. 

3. Classical Perturbations : 
3.1 Formalism 

There are three types of hnear cosmological perturbations: scalar, vector and tensor 
modes. The names refer to the way in which the modes transform under background space 
coordinate transformations (see e.g. Ref. 14). Tensor modes are gravitational waves, vector 
perturbations correspond to rotation and do not grow in time, and only the scalar modes 
couple (via the Einstein equations) to energy density and pressure. Hence, we shall restrict 
our attention to scalar type cosmological perturbations. 

The first step in the analysis of cosmological perturbations is to identify the gauge 
invariant combinations of Sgf^i,. The general scalar metric perturbation can be written in 
terms of four scalar functions (f), ip, B and E 

For simplicity, we have restricted our attention to the case of a spatially flat background. 
The following gauge transformations preserve the scalar character of 8 Qui,: 



(6) 



X* = X* + 7*-'^,^ 



where and ^ are functions of space and time. It is not hard to check that the induced 
changes of 0, -0, -B and E are 



a 

^ = ^ + 7^° (7) 
S = S + e°-^' 

E = E-^, 

where a prime denotes the derivative with respect to conformal time. Now it is a simple 
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exercise in linear algebra to find a basis of gauge invariant variables. A convenient choice is 

^ ^ (I) + a-^[{B - )a]' 



(8) 



a 



Note that in longitudinal gauge {B — E — 0) the gauge invariant variables become ^ — (j) 
and ^ = 

The second step of our analysis is to derive the equations of motion for the gauge invariant 
variables. In principle, this is straightforward. The linearized Einstein equations (2) are 
conveniently combined to yield equations for $ and ^. In practice, this computation is 
rather tedious unless a clever procedure is chosen. It is simplest^^ to consider first the 
transformation of the perturbation ^G^i^ of the Einstein tensor under (6), and to determine 
gauge invariant combinations (labelled with superscript (gi))'- 

SGf'^=SGU ^'^Go{B-E') 

5Gf'^ = 5G^ + - I Wg^) {B - E' ), (9) 

8Gf^'^ = 6G)+ ^^^G'J{B-E') 

where the background Einstein tensor elements are Evidently, the analogous combi- 

nations of 6Tlf are gauge invariant. Thus, the linearized Einstein equations can be written 
as 

^ g^Q gj^i^igi) _ (10) 

In this form, all the gauge dependence automatically drops out, and we obtain the following 
set of equations written exclusively in terms of gauge invariant variables: 



-m{H^ + ) + v^* = at^go^stI 

{n^ + -^')^i^A-KGa^5ti 
{2n' + 7^2)$ + + ^" + 27i^' + - V^D 5] - -f'^D^kj 



2 

= -47rGa2(5rj 



1 .... (11) 



where D — ^ — '^ and Ti — a' /a. 
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An alternative way to derive the above equations^) is to work in longitudinal gauge and 
at the end replace and ip by ^ and ^ respectively (and similarly for the matter variables). 

3.2 Applications 

As a first application of the classical theory of cosmological perturbations we shall con- 
sider the example of perfect fluid matter given by the energy-momentum tensor 

T| = (£+p)K>-p5^, (12) 

e and p being energy density and pressure respectively, and the four velocity vector of the 
fluid. In general, the pressure is a function of both e and entropy per baryon s, and hence 

Sp = clSs + tSs , (13) 

with Cs being the speed of sound. If r = 0, we have a pure adiabatic perturbation. 
The perturbation of is given by 

(5To° = Se 

5Tf ^{eQ+PQ)a-Hui (14) 
57] = -5p5i , 

where subscripts denote background quantities. Since (5Tj is diagonal, it follows immediately 
from the third equation in (11) that $ = ^ . This, in turn, leads to a significant simplification 
of the equations of motion for the gauge invariant variables. Prom the first equation in (11), 
we obtain 

V^* - m^' - 37^2$ = 47rG'a2fe(»^) . (15) 

This is a generalization of the Poisson equation to which it reduces in the Newtonian limit; 
and hence we call $ the relativistic potential. 

Equations (11) can be combined to yield the following second order equation of motion 
for $: 

+ 37i (1 + c^) - c^V^^ + [2n' + (1 + ^(?s)n^] $ = A-KGa\5s (16) 

Por adiabatic perturbations, the source term vanishes. On scales larger than the Hubble 
radius, the spatial gradients can be neglected. Under these conditions, equation (16) can be 
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recast as a "conservation law" 

C = (17) 
where the dot denotes the derivative with repect to physical time and 

C = ^ + $ 18 

^ 3 1+w ^ ' 

with w = p/e. The quantity ( was first introduced in Ref. 15 (see also Ref. 16). The 
above conservation law is easily applicable to many interesting issues. First, we note that if 
the equation of state is constant, then $ remains constant (the second solution of (17) is a 
decaying mode). However, during a phase transition w may change by a large factor. In this 
case, equation (17) implies that the relativistic potential $ will also change by a large factor. 
This is one of the key points in the computation of density perturbations from inflation ■'^'''^ 

To correctly describe fluctuations from inflation, we must consider a second apphcation 
of the classical theory of cosmological perturbations, namely a model with scalar field matter. 
The matter action is 



(19) 



semicolons denoting the covariant derivative. The induced energy- momentum tensor of the 
scalar field (p is 

= - I i - ^(y') } 5^ • (20) 

If we expand (p{x,t) about a homogeneous background field (po{t) 

<^(x,t) =<^o(t) + 5<^(x,t), (21) 



then the perturbation of Tg at the linearized level becomes 

5Tf = a-Vo5^,i (22) 
As in the case of a perfect fiuid, ST"^ is diagonal and hence $ = ^ 
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Inserting thios result and (22) into the general equations (11) and combining the resulting 
differential equations, we obtain the following second order equation for $ 



+ 2 7i - 



Since for a scalar field 



(n - ^) ^' - + 2(n' -^n^^ = o. (23) 



. 2 

1 + w = ^ , (24) 
e 

we can, like for perfect fluid matter, rewrite (23) as a "conservation law" identical to (17) 
and (18) when considering scales much larger than the Hubble radius. 

3.3 Fluctuations in Inflationary Cosmology 

To demonstrate how easy it is to apply the gauge invariant theory of cosmology pertur- 
bations, we shall consider the evolution of fluctuations in inflationary Universe models^^""*^^-*. 
We first note from (11) that on scales smaller than the Hubble radius 

— Kx)^(l)"- 

The calculation of density perturbations proceeds as follows: by evaluating (25) at the time 
ti{k) (see Fig. 1) when the wavelength under consideration leaves the Hubble radius, we 
determine the initial value of $,$(ij(A;)). By integrating (17) and using the fact that $ 
vanishes at both ti{k) and tf{k), the value of $ at the time tf{k) when the scale reenters 
the Hubble radius can be calculated with the result 

mm - I ^ a . (26) 

The coefficient a is 4/9 for in the radiation dominated phase and a — 2/h during matter 
domination. Using (25), the value of ^{tf) determines the late time value of the amplitude 
of the fractional density perturbation. 

In order to evaluate the amplitude of perturbations (26) in inflationary Universe models, 
a quantum analysis of the generation of fluctuations is required. However, already a rough 
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order of magnitude estimate yields interesting results. Prom (25) and taking the energy 
density perturbation to be given by the Hawking temperature (i.e. 5e ~ H^) we obtain 

^{k) < 1. (27) 

e 

Inflation is driven by a scalar fleld (p. During inflation, the equation of state is dominated 
by the potential energy density of (p. However, (/? is rolling and therefore 

. 2 

l + w{ti)^'^. (28) 

On dimensional grounds, LfP' ~ if^ and hence $(t/) ~ 1. We conclude that the change in the 
equation of state leads to a drastic amplification of the initial quantum fluctuations. Suc- 
cessful models of galaxy formation require $(t/) ~ 10^^. Thus, without careful adjustment 
of parameters, inflationary Universe models predict too large perturbations^^) . 
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4. Quantum Perturbations: 

4.1 Motivation 

The classical analysis of fluctuations in inflationary Universe models gives good insight 
into why initially tiny inhomogeneities are amplified by a large factor between when they 
are produced in the de Sitter phase and when they reenter the Hubble radius at late times. 
It is, however, only a quantum analysis which explains the origin of these perturbations. It 
is vacuum quantum fluctuations which are the source of the classical inhomogeneities which 
form the seeds for galaxy and cluster formation. 

A second motivation for considering the quantum theory of cosmological perturbations 
comes from the general problem of particle production in expanding background space- 
times. The usual^°^ treatment which is based on quantizing matter fields on an unperturbed 
cosmological background is inconsistent since matter fluctuations are intrinsically coupled 
to metric perturbations via the Einstein equations. Hence, we need to quantize metric and 
matter fluctuations in a unified way. 

In fact, the quantization of linear cosmological fluctuations is not more complicated than 
the well known quantization of matter fields in an external background: it is a straightfor- 
ward application of canonical quantization^^) . In synchronous gauge, the quantization of 
fluctuations was first discussed in Ref. 17. 

Since we only wish to quantize the physical degrees of freedom, it is advantageous to use 
the gauge-invariant formalism. Since this method reduces the number of degrees of freedom, 
it also leads to a substantial simplification of the analysis. 

The first step in deriving the quantum theory of cosmological perturbations^"*^'^^) is to 
determine the action for the fluctuations in terms of the gauge invariant variables. In general, 
it would be wrong to simply start from the classical equation of motion for perturbations and 
interpret it directly as an operator equation. This would lead to wrong canonical momenta 
and to a wrong normalization of the field operator^'^^^ . 

4.2 Formalism 

In the following, we shall briefly summarize the quantum theory of cosmological pertur- 
bations. For simplicity, only models with scalar field matter will be considered. For hydro- 
dynamical matter the analysis is similar^) . The formalism also applies to highter derivative 
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gravity theories^'^^^ . We will follow the method of Ref. 21 (see Ref. 4 for more details). 

The first and most involved step in quantizing cosmological perturbations is to write the 
action for fluctuations in terms of gauge invariant variables only. We start from the action 

J d^xV^R + Sm, (29) 

where Sm is the action for the scalar field (f. Next, we insert into (29) the expansion of g^i, 
and ip about a homogeneous background solution g^J and ipo 

(0) , X 

Quiy = QfiJ + dg^,,, ^^^^ 

ip = (po + d(p 

and expand the result in terms of powers of small quantities to find 

S = So + S2S (31) 

where 5*0 is the action of the background solution and S2S is quadratic in perturbation 
variables (the linear terms vanish because we are expanding about a solution of the equations 
of motion). We now use the constraint equations to simplify the action and drop total 
derivative terms. After a significant amount of algebra one obtains the following very simple 
form of ^25': 



(32) 

where v is a gauge invariant combination of matter and metric perturbations 



and 



a { + ^4 ) (33) 



The result (32) has the same form as the action of a simple scalar field with time depen- 
dent square mass —z"/z. Note that although the details of the reduction of the action are 
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somewhat involved, the final result is no surprise. We have seen in section 3 that for scalar 
field matter there is only one independent gauge invariant metric perturbation variable. Via 
the Einstein equations this variable is coupled to the gauge invariant matter fiuctuations. 
Thus, this is only one independent variable which expresses in a unified manner both matter 
and metric perturbations. 

From this point on, the quantization prescription is straightforward canonical quantiza- 
tion. Prom we can immediately write down the canonical momenta. After imposing the 
canonical commutation relations, we expand the operator v corresponding to the classical 
field V in terms of creation and annihilation operators and a^: 

V = I <fk [e'^-v*Ma-^ + e-'^-vMat] ■ (35) 

The mode functions V]f{r]) satisfy the equation 



v'^+(k'--]vk^O. (36) 



z 



Since (36) is a harmonic oscillator equation with time dependent mass, there will be 
quantum particle production^'^) . Modes of (36) which have positive frequency at some initial 
time to are no longer pure positive frequency at a later time ti > to. This leads to time 
dependence of expectation values of physical operators. For example, ii \ipo > is the vacuum 
state at time to, and A^fc(ti) = a'^ {ti)a^ (ti) is the number operator at time ti defined in 
terms of the operator coefficients of the positive frequency modes at time ti, then 

<MNk{tiMo>y^Q- (37) 

The final step is to compute the expectation values of the operators which determine 
the r.m.s. mass fluctuation. If 5M/M{k) is the r.m.s. mass perturbation inside a sphere of 
radius , then 

(^)\,).,3(^|y(,)W^,3li«*0PW*i_ (38) 

where in the final step we have replaced the classical perturbation by the expectation value 
of the quantum operator evaluated in the vacuum state IV'o > at the beginning of inflation. 
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This prescription for taking the quantum to classical transition has been discussed in Ref. 
5 and references therein. When evaluated at the time of Hubble radius crossing tf{k) (see 
Fig. 1), then using the relationship (25) between fe^^*) and $ one obtains 

{w) (kMk)) - < Mmi'li^o > . (39) 
In turn, the gauge invariant potential $ is related to the variable v by 

k^^ = -A7rG ^f-V ■ (40) 
rC \zJ 

Hence, the computation of the expectation value in (39) reduces to a straightforward evalu- 
ation of the expectation value of v^. 

Combining (39) and (40), we find 

(^) ^ k'WkUm^ (41) 

where Uk{yf) are proportional to the expansion coefficients of the operator $: 

^(x, V) = ^ / d'k [ul{v)e^-^al + «fc(r;)e-^-^a+] . (42) 

Equation (41) relates the r.m.s. mass perturbation resulting from quantum vacuum fluc- 
tuations to the solution Ukifj) of the classical equation of motion. At this point, we have 
established a consistent unified treatment of generation and evolution of cosmological per- 
turbations. 

Evaluating (41) for a model of chaotic inflation^^) with potential 

V{ip) = -if^ (43) 
n 

yields 



where / is the Planck length and is the characteristic wavenumber of the cosmic microwave 

M 



background. For n — 2 (and setting X — m^), a correct value of ^ for galaxy formation 
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requires m ~ 10 nipi, for n = 4 the requirement is A ~ 10 . As mentioned at the end of 
Section 3, inflation thus requires careful adjustment of parameters if it is to give the required 
value of density perturbations. 

This completes our brief survey of the quantum theory of cosmological perturbations. 
Although the formalism has been developed for scalar field matter (and in particular applied 
to quantum fluctuations in inflationary Universe models), it is much more general. Also for 
hydrodynamical matter and in higher derivative gravity theories, the action for perturbations 
can be reduced to a form like (32), and the quantization then proceeds as in the example 
discussed^) . 

5. Conclusions and Discussion : 

We have summarized the gauge invariant theory of classical and quantum cosmological 
perturbations. It allows a consistent unified treatment of the generation and evolution of 
linearized fluctuations in inflationary Universe models. 

In Section 2 we argued that a gauge invariant analysis of classical perturbations is phys- 
ically unambiguous and technically straightforward. It eliminates the gauge ambiguities 
associated with gauge-dependent approaches. The coordinate approach presented here is 
probably the most simple way of deriving the equations of motion for the gauge invariant 
gravitational potential. It is action based and hence allows standard canonical quantiza- 
tion. A gauge invariant analysis of the quantum theory implies that only physical degrees 
of freedom are quantized. 

The formalism presented here is practical and can easily be applied to problems of real 
cosmological interest. Already the classical theory leads to a useful "conservation law" (see 
(17) and (18)) which allows us to track the amplitude of perturbations on scales much larger 
than the Hubble radius in a very simple manner. 

The quantum theory of cosmological perturbations relates the expectation values of two- 
point functions which determine the r.m.s. mass fluctuations to mode functions which obey 
the classical equations of motion for the gauge invariant gravitational potential. This allows 
a unified analysis of the generation and evolution of density perturbations in inflationary 
Universe models (see (41)). A general formula for the amplitude of the resulting fluctuations 
is given. 
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In Ref. 4 we have performed detailed calculations of the spectrum of density perturba- 
tions in models with scalar field matter, hydrodynamical matter, and in higher derivative 
theories of gravity. The formalism also allows a discussion of entropy perturbations, it can 
be used to yield a simple proportionality between the microwave background temperature 
anisotropies and the gravitational potential and it can be applied to the generation and 
evolution of gravitational waves. 
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